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Abstract 

In this paper we ivestigate Stackel transforms between different classes of parameter-dependent 
Stackel separable systems of the same dimension. We show that the set of all Stackel systems of 
the same dimension splits to equivalence classes so that all members within the same class can be 
connected by a single Stackel transform. We also give an explicit formula relating solutions of two 
Stackel-related systems. These results show in particular that any two geodesic Stackel systems 
are Stackel equivalent in the sense that it is possible to transform one into another by a single 
Stackel transform. We also simplify proofs of some known statements about multiparameter Stackel 
transform. 
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1 Introduction 



Stackel transform is a functional transform that transforms a given Liouville integrable system into a 
new integrable system on the same Poisson manifold. It was first described by J. Hietarinta et al in [I] 
(where it was called the coupling-constant metamorphosis) and developed in 2 . It has been applied 
in ED U for the purpose of classification of superintegrable systems in conformally flat spaces. In 
[7J [5] the author described Stackel transform as a canonical transformation on an extended phase space. 
Applied to a Stackel separable system, this transformation yields a new Stackel separable system, which 
explains its name. 

Originally, only one coupling constant, entering linearly in one of the Hamiltonians of the system, was 
used. In paper [9] a multiparameter generalization of Stackel transform has been introduced. This gener- 
alization allows for a nonlinear dependence of Hamiltonians of the system on several coupling parameters, 
thus much enlarging the class of admissible Stackel transforms. Also, this generalized transform results 
in a class of reciprocal transformations that has been applied in [TTJ for analyzing weakly-nonlinear 
scmi-Hamiltonian hydrodynamic-type systems. This indicates that Stackel transform is a useful tool 
for studying various integrable systems. It can also be generalized for studying systems of ODE's of 
evolutionary type with integrals of motion, see [11] . 

In this paper we use the approach developed in and further in [T^] to show that all Stackel systems 
of the same dimension n can be split into equivalence classes such that every two members of the same 
class are Stackel equivalent in the sense that there always exists a single n-parameter Stackel transform 
between arbitrary two such systems. In order to do this we consider Stackel transforms inside given classes 



1 



of Stackel systems, a problem not considered in previous papers. We also give an explicit, compact form 
of this transform, making the formulas more transparent then these in [5] and |12j . Also, we present a 
corresponding reciprocal transform between solutions of these Stackel-related systems confined to proper 
submanifolds of the phase space. We also clarify and in some cases also repair a number of formulas and 
simplify proofs of a number of statements given in [9] . Two extensive examples are given at the end of 
the paper. 



2 General Stackel transform 

In this section we present some facts about multiparameter Stackel transform. Consider a manifold M 
equipped with a Poisson tensor II. Denote the space of all smooth functions on M by C°°(M). The 
mapping {•,•} : C°°(M) x C°°(M) -> C°°(Af) given by {f,g} u = (df,ILdg) (where (•, •) is the dual 
map between cotangent and tangent spaces) is called Poisson bracket and it turns C°°(M) into a Lie 
algebra. Suppose we have r functions (later: Hamiltonians) hi : M —¥ R on M, each depending on k < r 
parameters ai, . . . , a k so that 

hi = hi(x,ai,...,ak), i = l,...,r, (1) 

where x G M. Let us now from r functions in ([1]) choose k functions h Si , i = 1, . . . , fc, where {si, . . . , Sk} 
C {1, . . . , r}. Assume also that 

det(dh s Jda 3 ) ^0. (2) 
The condition ^ means that the system of equations 

h Si (x,ai,...,ak) = &i, i = l,...,k, (3) 

(where <Sj is another set of k free parameters, or values of Hamiltonians h Si ) involving the functions h Si 
can be locally solved for the parameters a, yielding 

CH = h Si (x, a k ), i = l,...,k, (4) 

where the right hand sides of these solutions define k new functions h Si on M , each depending on k 
parameters c^. Finally, let us define r — k functions hi with i — 1, . . . , r and such that i ^ {s\, . . . , s&} 
by - in accordance with Q - substituting h Si instead of ai in hi for i ^ {s\, . . . ,Sfc}: 

^ = ft ila 1 ->h. 1 ,...,a»->hV i = 1 '-"' r ' it {81,. -.,8k}. ( 5 ) 

Definition 1 The functions hi — hi{x, a%, . . . , &k)t * = 1 , - - - , defined through |^) and |3P are called 
the (generalized) Stackel transform of the functions {TJ) with respect to the indices {s\, . . . , Sk} (or with 
respect to the functions h Sl , . . . h Sk ). 

Note that unless we extend the manifold M this operation can in general not be obtained by any 
coordinate change of variables. It is also easy to see that if we perform again the Stackel transform on 
the functions hi with respect to h Si we will receive back the functions hi in (p}. In this sense the Stackel 
transform is a reciprocal transform. Note also that neither k nor r are related to the dimension of the 
manifold M. 

Example 2 The simplest situation occurs when k = r = 1. Consider, after fl^, the Fokas- Lager strom 
potential on the four- dimensional phase space M with coordinates (x,y,p x ,p y ) : 

h=\{ P l+ P 2 y )-\a{xy)-^ 

Solving the equation h = a with respect to the only parameter a (called in J7]/ a coupling constant) one 
obtains 

a=\{xyf"[pl+pl) -\{xyYl*a = h 
which can be shown [1] to be equivalent to the axially symmetric potential p A . 
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Stackel transform has two important properties that make it well suited for study of integrable systems: 
as we will see in Theorem [31 it preserves functional independence and it also preserves involutivity with 
respect to the Poisson tensor II. Moreover, as it will also be demonstrated in this paper, it maps a Stackel 
separable system into a new Stackel separable system which explains the name of this transformation. 

In the special but nonetheless important for this paper case when functions (JlJ depend linearly on 
parameters ai it is possible to write down the Stackel transform explicitly. Suppose therefore for the 
moment that the functions in ([!]) have the form 



r C0 



» = 1, 



(6) 



hi = Hi + ®j H ; 

J'=l 

The equations (J3J defining the first part of the Stackel transform take then the form of a system of k 
linear equations in k unknowns a\ 



with the Cramer solution for a,; = h s 



where 



, Oik 



ai 



3 = 1 

of the form: 

h St = det Wi/det W, 



(7) 



W = 



H.. 



(i) 



H.. 



(i) 



H 



H. 



(fc) 



fe) 



is the A; x A; matrix det {dh s J dctj) given in ([2]) (so that det W ^ 0) and where Wi are obtained from W by 

(i) 

replacing Hs, in the i-th column by &j — H Sj for all j — 1, . . . , fc. The second part of the transformation, 
i.e. formulas ([5]), reads now 

k 

hi = Hi + Y]h s H^\ i = l,...,r, i <£ {si, . . . , s k } 



where h Si are given by For k = 1 the above transformation reproduces the original Stackel transform 
presented in [T] and [2]. 



3 Stackel transform for completely integrable systems 

Let us now discuss the Stackel transform and the corresponding reciprocal transform between two Liouville 
integrable systems. Suppose therefore that dimM = 2n and that we have exactly n (so that r = n now) 
functionally independent functions (Hamiltonians) 

hi = hi(x, ai, ■ ■ ■ , oik), i = l,...,n 

that depend on k < n parameters ctj and that are for all values of ai in involution with respect to 
a nondegenerate Poisson bracket II: {/i;,/ij} n = for all These functions yield n commuting 

Hamiltonian systems on M: 

dx 

— =Hdhi=Xi,i = l,...,n (8) 

dU 

so that Xi are n commuting Hamiltonian vector fields on M. Consider now a new set of n functions 
(Hamiltonians) hi obtained from hi by a Stackel transform performed with respect to h Sl , . . . , h Sk (recall 
that k < n). These functions define a set of Hamiltonian flows on M, the vector fields of which are given 

by 

dx ~ — • 

-^r = Udhi = X,^ i = l,...,n. (9) 
dU 
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Let us now consider the relation between the Hamiltonian systems ([5]) and ©■ As it follows from ^ 
and ([3]) the following identity is valid on M and for all values of parameters af. 

h Si (x, h Sl (x,a 1 ,...,a n ),...,h Sk (x,ai,...,a n )) = a l , i = l,...,k (10) 

Moreover, the second part of the transformation, i.e. formula ([5]) can be written as the following identity 
on M (valid again for all values of Si): 

hi(x,ai,..., S„) = hi(x,h si (x, Si, . . . , S„), . . . , h Sh (x, Si, . . . ,S n )), i = 1, . . . ,n, i ^ {si, . . . , s/J 

(11) 

Differentiating (fTO)) with respect to x p we find that 

k dh ~ 

'II'- Y,-^-'"'-- i = l,...,* (12) 

2=1 J 



while differentiation of (jTTJ) gives 



A dh 



dh i = dhi -^2g^-dh Sj , i = l,...,n, i £ {s u ■ ■ ■ ,s k } . (13) 

2=1 J 

The transformation (|12p ~ (|13p can be written in a matrix form as 

dh = Adh (14) 
where we denote dh = (dhi, . . . , dh n ) T and dh = (dhi, . . . , dh n ) T and where the n x n matrix A is given 

by 

Aij = Sij for i £ {si, ...,8k}, j $. {si,...,s fe } 

A, 4 j = for i = 1, . . . , k, j $ {s%, ...,s k } 
A SiSj = - d -^io,i,j = l,...,k. 

OCtj 

^From the structure of the matrix A it follows that 

'dh s 



det A = ± det 



da 3 



so that det A ^ due to the assumption (J5J). Thus, the relation P^)l can be inverted yielding dh — A 1 dh. 
This leads to an important theorem [5], mentioned in Section 2. 

Theorem 3 1. If the functions hi are functionally independent then so are hi. 2. If the functions hi are 
in involution with respect to the Poisson tensor II (for all values of at), then the functions hi are also in 
involution with respect to II for all values of oti . 

Proof. 1. Assume that hi are functionally independent. It means that dhi are linearly independent 
at every point of M. Then, by (fT4"| and by the fact that det A ^ the differentials dhi also are 
linearly independent and thus hi are functionally independent on M. 2. Assume {hi,hj} u = for all 
i,j = l,...,n. Then 

(n n \ 

n n 

= E (A- 1 ) tfi (A- 1 ) JIa (dft, 1> ndft, a )= E (^M^W^-^in-o 



h,h = l Ji,/a=l 
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Theorem [3] implies that the system © is again Liouville integrable so that Stackel transform maps 
a Liouville integrable system into a Liouville integrable system and from its proof it also follows that 
Stackel transform also preserves superintegrability [3]- [6]. 

Since — Hdhi and Xi = Hdhi we obtain from (TT2)) - ([T5)) that the Hamiltonian vector fields X t = 
Hdhi and X, = Tldhi are related by the following transformation 




This means that the hamiltonian vector fields Xi and Xi span the same n-dimensional distribution 
on M and also that the vector fields X Si and X Si span the same fc-dimensional subdistribution of the 
above distribution. The transformation (|15 p -(|16 [) can be written in matrix form as 

X = AX (17) 

where we denote X = (Xi, . . . , X n ) T and X = (X\, . . . , X n ) T and where the n X n matrix A is given 
above. 

All the vector fields Xi and Xi are naturally tangent to the common level surface M a ^ of all the 
Hamiltonians h Si : 

JW ,5 = {x £ M : h Si (x,ai,. . . ,a k ) = on, % = 1, . . . k} 

(since they are obviously tangent to common level surfaces of all n Hamiltonians hi) so that if Xq £ M a & 
then the multiparameter (simultaneous) solution 

x = x(tx, ■ ■ ■ ,t n ,x ) (18) 

of all equations in JSJ starting at xq for t = 0, will always remain in M a ^ and the same is also true for 
multiparameter solutions of ©• Note that the surface M a 5 depends on the choice of 2k parameters 
and Sj and that its codimension is k so that dimA/^ 5 = 2n — k > n. Note also that, by virtue of the 
definition of the Stackel transform, M a & — M 5 a where 

M S:a = (x G M : h Si (x,ai,. . . ,a k ) = on, i = 1, . . . k\ 

The relations p3|) - (TT6]) can be reformulated in the dual language, that of reciprocal (multi-time) 
transformations. 

Theorem 4 The reciprocal transformation ti = U(ti, . . . , t n , x), i = 1, . . . , n given by 

dt = A T dt (19) 

(where dt = (dt\, . . . , dt n ) T and dt = (dti, . . . , dt n ) T ) transforms the n-parameter solutions i f_?<S|) of the 
system (0i to the n-parameter solutions x = x(t\, . . . , t n , Xq) of the system (0i (with the same initial 
condition x(0) — xq) in the sense that for any xq G M we have 

x(ti(h, . . . ,t n ,x ), ■ ■ - ,t n (ti, ...,t n , x ),x ) = x(ti,. . . ,t n ,x ) 

for all values of ti sufficiently close to zero. 

The transformation (|19p is well defined since the right hand side of (TH)|) is an exact differential, as it 
follows from the above construction. It means that it is possible (at least locally) to integrate (TT9|) and 
obtain an explicit transformation ti — U(ti, . . . , t n , x) that takes multi-time (simultaneous) solutions of 
all hamiltonian systems ^ to multi-time solutions of all the systems in ^ . 
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In a specific but important for us case when k — n (i.e. when the number of parameters and the 
number of hamiltonians coincide so that the Stackel transform consist only of the first part i.e. @), the 
matrix A simplifies insofar as 

a dhi . . 

so that the formulas (IT5|) - (IT51) simplify to the single formula 



j=i j 



, n, 



while (|19[) can be explicitly written as 

" dh 

dii = -^2 Tr^dtj, i = 1, . . . , n. (20) 
i=i aai 

and our manifolds M a 5 become in this case level surfaces for all the hamiltonians hi(x, a) and also level 
surfaces for all the hamiltonians hi{x,a). 



4 Classical Stackel systems 

Consider a set of Darboux coordinates (A, fi) = (Ai . . . , , A n , fii, . . . , /ii) on our 2n-dimensional Poisson 
manifold M equipped with a Poisson operator II (so that II = X)i<j ax" ^ af - )- A classical Stackel system 
is a system of n Hamiltonians Hi on M (that originally do not depend on any additional parameters a 
so that they can not be a subject of any Stackel transformation) originating from a set of n separation 
relations [T3] of the form: 

n 

cx(Ai)+^W =/(A i )/i?, i = l,...,n, (21) 
j=i 

where / and a are arbitrary functions of one argument and where all ji € Z, i = 1, . . . , n, and are such 
that no two 7^ coincide. Thus, a particular Stackel system is defined by the choice of integers 71, . . . , 7„ 
and by the choice of functions / and a. Customary one can also treat this system of relations as n points 
on (n copies of) the following separation curve 

n 

P(X,H) ee .7(A) + 5>,A^ = ,/(A)/i 2 , (22) 
3=1 

in \[i plane which helps us to avoid writing too many indices. The relations (f2"Tj) (or n copies of (|2"2"10 
constitute a system of n equations linear in the unknowns Hi. Solving these relations with respect to 
Hi we obtain n commuting (since the right-hand sides of formulas (|2ip commute) with respect to II 
Hamiltonians (known in literature as Stackel Hamiltonians) on M of the form 

Hi = 1^ KiG[i + Vi(X) i = l,...,n, (23) 

where we denote A = (Ai, . . . , A„) T and /i = (fix, . . . , fi n ) T ■ The functions Hi can be interpreted as n 
quadratic in momenta /j, hamiltonians on the phase space M — T* Q cotangent to a Riemannian manifold 
Q (so that Ai, . . . , A„ are coordinates on Q) equipped with the contravariant metric tensor G depending 
on function / and the choice of the set 7. The objects Ki in (1231) can be interpreted as (1, l)-type 
Killing tensors on Q for the metric G. The metric tensor G and all the Killing tensors Ki are diagonal 
in A- variables. Note that by the very construction of Hi the variables (A, //) are separation variables for 
all the hamiltonians in (|23|) in the sense that the Hamilton- Jacobi equations associated with all Hi = ai 
admit additively separable solutions W — Y17=i^i(^i' °)- 
The relations (|2ip can be written in a matrix form as 

S~H = U 



G 



where H = (Hi, . . . , H n ) T , and where U is a Stackel vector of the form 

U = (/(Ai) M f - tf(Ai), . . . , /(A„)^ _ a(A„)) T , (24) 



while the matrix is a classical Stackel matrix of the form 



Note that our assumption that no 7; coincide means that det(S' 7 ) 7^ 0. Thus, the hamiltonians 
be obtained in a matrix form as 

H = S^U, 

which also means that the metric G in (|23l) can be expressed as 



G = diag (/(AO (5 7 - 1 ) n , . . . ,/(A„) (S~ l ) ln ) , 

so that the Killing tensors Ki in (|2"3"| arc 

Ki = diag ((S; 1 ^ I (S-') 11 , . . . , (S^) m I (S^) ln ) ,< = !,...,» 

(note that K\ = I). Let us now turn our attention to the scalar functions Vi : Q — >• R in (1231) . Relations 
(|221) and (|23|) imply that V^(A) satisfy the following separation curve 

(7(A) + V1A 71 + V 2 X~< 2 +... + KA 7 " = 0, (26) 

so that they depend on the choice of integers 7, and the choice of the function a. We will therefore 
denote them as V, . In case when a(X) is a monomial, i.e. when a(X) = X k with k € Z, Vi depend on 
k and they will be denoted by (so that — V^ X ■*) to shorten the notation. Thus, the potentials 
v} k \X) (they still depend on all 7,) satisfy the following separation curve 

A fc + V-^A 71 + y 2 (fe) A 72 + . . . + V^X^ = 0, (27) 
which in matrix form can be written as 

S 7 V (fc) = -A fe (l,.-.,1) T , (28) 
where = (v} k \ . . . , V^ k) ) T and A =diag(A x , . . . , A„). This means that 

VM = -S;\1,...,1) T 

so that = S^AS-yV^, = S^A 2 S 7 V^ = (S^AS^S^AS^V^ and so on. Similar argu- 
ment applies also for negative k. Thus, denoting 

F 7 = 5" X AS 7 (29) 
we get the compact formula for the potentials (presented first in |12|): 

It is now an immediate consequence of the above formulas that for any meromorphic function a(X) we 
have 

— a(F 1 )V^°\ (31) 

where = (v} a \. . . , V^ ] ). The matrix F 7 given in (|2"91 has been called control matrix in [T3] where 
it appeared in the context of quasi-bi-Hamiltonian representation. Notice also that if the system (|2"21 
is normed by 7„ = 1, then the potential attains a particularly simple form V' ) = (0, ... ,0, — 1) T . 
This follows immediately from (|28l) . 

Now, by writing down the inverse Jacobi problem for all the Hamiltonians (f23|) we can arrive at the 
following remark that will be useful in the next section when we discuss reciprocal transforms between 
different Stackel systems. 
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Remark 5 On the level surface M a = {i£ M : Hi = G R} the mutliparameter (multi-time) solutions 
\i = Aj(ii, . . . ,t n ,xo) of all Hamiltonian systems defined by the separation curve h22\) or equivalently by 
all Hamiltonians \23)) attain the following Abel-Jacobi differential form 



dt = s r E - dX =, 

' y/f(X)P(X, a y 



(32) 



where dX/ \J /(A)P(A, a) means a column vector with components dXi/ \J /(Ai)P(Aj, a). 



Note that solutions (j32j) define in a standard (canonical) way the corresponding multi-time solutions 
for the momenta p,i — fJ-i{ti t . . . , t n , xq). 

A particular subclass of Stackel systems is given by choosing % = n — i. The separation curve (|2"21 
attains then the form 

n 

a(A)+5>,A-->=/(A)M 2 
j'=i 

and the originating hamiltonians Hi constitute a completely integrable system that we call a system of 
Benenti type due to S. Benenti's contribution to the study of these objects [T3],[I5]. In this case it is 
possible to give compact formulas for many objects introduced above. Thus, for example, the metric G 
and the Killing tensors Ki in (1231) are given explicitly as 



G = diag 



( f(M 



V Ax 



/(An) 



r- , I d Pi <>/>. 

K l = - diag — — ,• • • , — - 



Ai = n (A* 



1. 



A,) 



Here and below pi 



K dXi' , dX ri/ 

Pi (A) are Viete polynomials (signed symmetric polynomials) in A: 



MA) = (-i) 



A, 



l<Si <S2<-..<Si<n 



A s 



i = 1, . . . , n 



(33) 



that can also be considered as new coordinates on the Riemannian manifold Q (we will then refer to them 
as Viete coordinates). Notice again that the Killing tensors Ki do not depend on a particular choice of 
/ and a. It can be shown that as long as / is a polynomial of degree < n then the metric G is flat while 
if / is a polynomial of degree n + 1 then G has constant but non-zero curvature. For the Benenti class 
the control matrix F 7 in (|29[) attains the simple form 



F = 



-Pi 

-P2 



V -Pn 



1 

o / 



(34) 



and since V( Q ) = (0, 0, . . . , 0, -1) T we easily obtain that the potentials = FV^ =(0,0,... 0, -1, 0) T , 
y{2) = F 2y(0) = (o,0,...0,-l,0,0) T up to V^' 1 ^ = F n ~ 1 V^ = (-1,0, ...,0) T , are trivial (con- 
stant), = F n V(°*> = (pi,...,p n ) is the first nontrivial positive potential while = F~ X V^ 
— (l/p n , pi/p n , . . . , p n -i/p n ) T and so on. More information on Benenti systems can be found in 



5 Stackel equivalence of Stackel systems 

We will now turn to the main question of this article: how to relate two Stackel systems by a single Stackel 
transform and in such a way that their solutions are related by a reciprocal transform? As we mentioned 
above, The Hamiltonians Hi defined by (|2Tj) or by (|22]l do not depend on any additional parameters 
ai so in order to perform a Stackel transform on (|2Tj) we have to embed it into a parameter-dependent 
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system. Of course, there is infinitely many ways of embedding of our Stackel system into an n-parameter 
system but the choice below is natural in the sense that the corresponding Stackel transform transforms 
a Stackel system into a new Stackel system. Thus, consider n Hamiltonians hi — hi(X,fj,,a) defined by 
the separation curve 

n n 

P(X, h, a) = <j(A) + ^hjX 1 ' +R~ 1 Wj2 a i xSl =f( X )^ ( 35 ) 
i=i i=i 

where 71 , . . . , 7„ and Si,...,S n are two sequences of integers such that no two % coincide and similarly 
no two Si coincide (but we do admit the possibility that some or all of 7, coincide with some Si) and 
where -R(A) is an arbitrary meromorphic function of one variable so that 

R(X) = l[(X-l3 s )f[(\-l3' s r 1 

s=l s=l 

for some (complex in general) constants fi\ , . . . , (3^ and /3[ , . . . , /3' k . This function can be generalized to 
a matrix function i.e. we define, for any n x n matrix A (A-dependent or not) 

R(A) = Y[(A-f3 s )l[(A~f3i)- 1 (36) 



(note that all the terms in (|5r?)) commute so that there is no ordering problem here). The relations (|55|) 
can now be written in a matrix form as 

S-yh + R- 1 (A)S s a = U (37) 

where S* 7 and S$ are two Stackel matrices given by 425]) (so that (5 7 ) ij = XJ* and = A* 3 ), 

/i = (hi, . . . , /i n ) T is the column vector consisting of Hamiltonians hi, a = (ax, . . . , a n ) T , U is the column 
vector given in and where A = diag(Ai, . . . , A„) as before. Solving (|3"T1) with respect to h we obtain 

h = S~ X V - S- 1 R- 1 (A)S 5 a = S~ X U - S- 1 R^\A)S 7 S^ 1 S s a. (38) 

Lemma 6 In the notation as above 

S- 1 R(A)S. y = R(Ey). 

Proof. We show it for R(X) = X — /3 as the general statement follows easily by developing the argument 
below. 

S' 1 (A - p) S 1 = 5- 1 A5 7 - 01 = F 1 - )3I = R(F 1 ). 

■ 

Thus, introducing the shorthand notation 

Ws^ — — S y 1 Sg 

we see that (|38|) can be written as 

h = H + R- 1 (F 1 )W Sn a (39) 

where H = S~ 1 U is the part of h that is independent of parameters an (cf. ((6])). Let us shortly analyze 
the structure of the matrix W& n = —S~ 1 Ss- Assume that cr(A) in (|26|) is a polynomial of the form 
a W = S"=i £i^ Si - Then, as it follows from (f!M| and from the definition of potentials V^ k ' : 

n 

1=1 

On the other hand, the formula (|26p can now be written as 

S^V^ +S S £ = 
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so that V (a) = -S^SgZ = W s ,~£ which implies that 

where V^O — i^ j y(°) in accordance with (|30l) . Therefore, the formula ([39]) can be written as 

n 

hi =Hi+J2 (fl -1 ^))^ * = l,...,n. (40) 

Let us now perform an n-parameter Stackel transform of the system given by the curve Q35p . Since 
the number of parameters oti and the number of Hamiltonians hi are both the same (= n) the Stackel 
transform consists only of part (j4]) and is therefore generated by the relation h — a (which implies h = a) 
in the vector notation as above. We are now in position to formulate the main theorem of this paper. 

Theorem 7 The n-parameter Stackel transform generated by h — a transforms the set of n Hamiltonians 
h defined by 135\) into the following set of Hamiltonians 

h = -W^R[FJH + W^R(F y )a (41) 

(where h = (hi, . . . , h n ) T and similarly a — (5i, . . . , a n ) T ) which constitute a new Stackel system with 
the separation curve of the form 

n n 

P(X,h,a) = R(X)a(X) + J R(A)^5 j A^ + ^ 3 A^ = R(X)f(X) f i 2 . (42) 

Moreover, the reciprocal transformation 

dt = -R~ 1 (Fg)Wg t7 dt = -W^R-^F^dt (43) 

transforms n-time solutions x = x(t\, . . . ,t n , Xo) of the system \35)) into n-time solutions x = x(t\, . . . ,t n , xq) 
of the system J^£| ). 

Note that in spite of the fact that we introduced both systems in the (A, /z)-variables the matrix 
formulas (|4Tj) and (|43|) are not tensor and that they are coordinate-free. They can be therefore freely 
applied in any coordinate system on M , which will be used in the examples further on. 
Proof. Multiplying the curve ff35|) by R(X) we obtain 



R(X)a{X) + R(X)Y,hjX^ + = R{X)f(X)n 2 

i=i o=i 

which after the Stackel transform h = a (so that h = a) obviously attains the form (|42"T) . Let us therefore 
show the formula (|41[) . The separation relations implied by (|42[) can be written in matrix form as 

R(A)Sja + S s h = R(A)U 

with the column vector U as in (|24j) . Solving this with respect to h we obtain 

h = S^R^U - Sg 1 R(A)S J a = Sg 1 R(A)S J H - Sg 1 R(A)S 1 a 
= (S^R(A)Ss) (S^Ssy 1 H - (Ss 1 R(A)S s ) (S^Ss)' 1 a 
= -R(F S )W^H - R(F s )W^a 

so the only remaining thing is to show that R(Fs)Wf* = W^} j R{F 1 ) which is equivalent to the statement 

W Sn R(F s ) = R(F 7 )W 5n 
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that can easily be proved in a fashion similar to proof of Lemma |51 Finally, the formula (|4"3")l is obtained 
by inserting (|35J into 



dt = ~\da) dt = ~ \ R (^Ww dt ' 

where we use the fact that R{A) T — R(A T ). m 

Let us also remark that the relations (j4Tj) can be explicitly written as (cf. (|40|l ): 

n n 

hi = - R ( F skVj H k + J2 R(FshV^ k) a kl i = 1, . . . ,n 

j,k—l j,k—l 

where the potentials Vj are defined by the separation curve (cf. (|27|l ) 

X k + V} k) X Sl + V 2 (k) X 52 +... + X Sl = (44) 

so that 

y(k) = F ky(0)_ 

Notice that due to the form of and the potentials and are related by the Stackel 
transform — a, — a. 

Let us also present an alternative way of proving the formula (|43l) , directly involving solutions of (1351) 
and It follows from Remark [5] and from the above considerations that the multi-time solutions of 

the systems (|3"5|) and (14"21) on any common level surface M a s attain the form 

dt = ST - dX = , dt = Sj , dX (45) 



y/f(X)P(X,a,a)' ^R(X)f(X)P(X,a,a) 



Now, it is easy to see that P(X,a,a) — R(X)P(X,a,a) so that, by ([45} and by the fact that R(A) is 
symmetric 

£= S J , dX = Sj ■ dX = SjR-\A) - dX = 

\jR(X)f(X)P(X,a,a) v^ 2 (A)/(A)P(A, 5, a) V/(W«» 

= SfR-\A) (S^dt = (S- 1 R-\A)S s f dt = (S^SsS^RrH^Ssfdt = -(W^ir^f* 
thus yielding 

«£= -R-^F^W^dt = -W^R-^F^dt, 

which is what we wanted to prove. 

On the level of Stackel transforms Theorem [7] leads to the following corollary: 

Corollary 8 Assume that fi(X) ^ and /2(A) ^ 0. Any two Stackel systems of the form 

n 

a 1 (X)+J2HjX^ = /i(A) M 2 

j'=i 

n 

a 2 (X) + ^2H j X 5 ' =/ 2 (A) M 2 
i=i 

£/ia£ satisfy the condition 

/ 2 (A)a 1 (A)=<T 2 (A)/ 1 (A) (46) 
are Stackel-related by the single Stackel transform 

H = -W s -iR(F 7 )H (47) 

im'f/i i?(A) = ^ 7^ . In particular, any two geodesic Stackel systems (i.e. with o\ = a 2 = 0) are connected 
by the Stackel transform f^9j ) with R(X) — ^ ^ . 
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Proof. By Theorem [TJ the Stackel transform (l47t with R(X) = f^x] transforms the first of the above 
Stackel systems into the second one provided that -j^y = ^l^r which is exactly the condition P5|) . In 
case of geodesic systems, the condition (j46]l is always satisfied. ■ 

Proposition 9 The condition \^b^ splits all Stackel systems of the form \21\) into equivalence classes 
since it is an equivalence relation. 

Proof. Indeed, if £g} = = i?i(A) and g$ = f$ = i? 2 (A) then = f$ = i? 2 (A)i? 1 (A) so 

this relation is transitive. Further, if ^jjy = ffjiy = then jjjjy = fjjxy = "r^aJ so ^ n ^ s relation is 

reflexive. Finally, ^jj) = §7(iy = 1 so it is a symmetric relation. 
■ 

Our formulas contain two special cases: when 7 — 5 and when R = 1. In the first case (i.e. when 
7i = 5j, i = 1, . . . ,n; Benenti systems are in this class) we relate systems belonging to the same class, 
where the class is understanding as a fixed sequence 71, j n , and differ by / and c.The matrix W 7l7 = — I 
while F$ — F 1 so that the formula (|41"]) becomes 

h = R(F 1 )H - R(F 7 )a (48) 

while (j43|) attains the form 

dt = R- l (F*)dt. 

In the second case (R — 1) we relate systems from different classes, i.e. (71,..., 7„) and (5i,...,S n ) 
respectively, which share the same / and a. The formula (|41|) becomes 

h=-W£H + W£a 

while the formula (|43|) attains the form 

dt= -Wj^dt. (49) 

Thus, the general transformation between the systems (|55|) and (|4"2"]l can be considered as composition 
of two transformations: a map between two Stackel systems from the same class (i.e. with 7 = S) but 
with different / (i.e. metrics) and/or different a and the transformation between two Stackel systems 
sharing the same / and a but from different classes. Both these transformations commute. 



6 Examples 

We will now present two examples of our formulas. In order to relate to known integrable systems we 
will present both examples in their natural (physical) coordinates. 

In our fist example we will relate two families of separation curves for n = 2, namely 

P(X, h, a) = (t(A) + M + h 2 + X(a 1 X 2 + a 2 ) = ^A/i 2 (50) 

and 

P(X,h,a) = (j(A)A- 1 + A~ 1 (5iA + 5 2 ) + hiX 2 + h 2 = ^ (51) 

which are particular cases of (|35|) respectively (|42|) with (71,72) = (1,0), (Si, S%) = (2,0) and with 
R = A -1 , while cr(A) is for now assumed to be an arbitrary rational function of A. The family (f50| 
contains in particular a well known Henon-Heiles (HH) system while the family (fBTj) contains in particular 
Drach system {2D]. Note that in both of the above curves one 7$ coincides with one Si and therefore we 
can regroup the terms in the above curves to obtain 

ct(A) +aiA 3 + (hi +a 2 )X + h 2 = ^A^ 2 (52) 

for the HH family ([50]) and 

a{X)X- 1 + hiX 2 + (h 2 + Si) + 5 2 A- X = i^ 2 (53) 
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for the Drach family (jSTj) . Since both metrics are fiat we will now construct both systems in their 
respective flat coordinates. As an intermediate step, we will write the systems in Viete coordinates (1331) 
that now attain the form 

Pi = -Ai-A 2 , /52 = AiA 2 . (54) 
In the above coordinates the control matrix F 7 (cf. (|29|) and ([34"]) ) of the HH family is 



F 1 = 



-Pi 1 
-pa 



and as we remember it is coordinate free. The passage to the flat coordinates (xx,x 2 ) is given by the 
point transformation [21] 

1 2 

Pi = xx, pi = —^x 2 . 
The metric G and the Killing tensor K 2 are now 



1 \ T ^ ( -\x 2 



so that geodesic (i.e. with cr(A) = 0) Hamiltonians of ({52} are 

1 2 1 2 1 2 1 

= 2^1 + 2^ 2 ' 2 = 2 :El2/2 ~ 2 X2VlV2 ' 

where (j/i, j/2) T are momenta conjugate to (xx,x 2 ). Note that (xi,x 2 ) are not only flat but also orthogonal 
coordinates for (|52|) . The matrix F 1 is now 



— X\ 1 

1 H o 



so that the Hamiltonians of (f5"2j) in the flat coordinates are 

hx = Ex + v}^ - (x 2 + ax - a 2 

hi = Ex + + -xia^ai, 

where due to (|3"Tj) we have VW = a(F y )VW = ct(F 7 )(0, -1) T . Let us now confine ourselves to the 
following three-parameter set of separable potentials: 

ct(A) = M 5 + b 2 X 4 + b 3 X 2 , b t G R. (55) 

We can now calculate V*"' explicitly and we find 

hx=Ex- bx (x\ + ^x\x\ + + h (xf + -Xxxj) + b z x x - (^x\ + ^x?^j ax - a 2 

h 2 = E 2 + bx (yjx\x\ + ^xxx{) - b 2 {^x\x\ + ~ ^b 3 xj + ^xxxjax 

and for ax = a 2 = 0, bx = 63 = 0, b 2 = 1 we receive the classical Henon-Heiles system. 

Consider now the Drach family (|53|). The matrix F$ in Viete coordinates ([54)) has the form 



F S = 



pi pi 

£2. P2 

Pi Pi 



We pass now to flat but non-orthogonal coordinates (x, y ) given by [20] 

pi = -2x*, p2 = x -y. 
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The metric G and the Killing tensor K 2 of this system are 



\ \ ~ _ ( -(x + y) 2x 



G il oj'^^ 2y -(x + y) 
so that the geodesic Hamiltonians of ([S"3"]) are 

1 ~ 1 o 1 o 1, 

#i = -jPxPy, E 2 = -xp x + -yp y - - (x + y)p x p y , 
where (p x ,Py) are momenta conjugate to (x,y). The matrix F$ becomes 

\x~^ (3a; + y) ' 



Fs = 
and hence 



_ . „ , , , , - T 2 

2* V"* 1 ^ 2~ 

- y) 2 ^x~^(x-y) 



hx = E x + V} a) + ~x-i(x - y)- 1 a 2 

fr 2 = £ 2 + V" 2 (<T) - Si - - y) _1 (3a; + y)5 2 

where due to ([31]) and to ([53]) we have = <t(F s )F^ 1 V < - 0) = a(F s )Ff 1 (0, -Vf . For our particular 
choice of u(A) as in ([55]) we have cr(A)A _1 = &i A 4 + 6 2 A 3 + 63A and the Hamiltonians hi attain the explicit 
form 

hi =E t - 2b 1 (x + y)- ^b 2 x~ ^ (3x + y) - ^b 3 x'^ + ^x~i(x - y)^ 1 a 2 

h 2 = E 2 + h(x - y) 2 + -b 2 X~^(x - y) 2 - ^b 3 x~^(x - y) - Si - -x~^(x - y) _1 (3a; + y)a 2 . 

The above Hamiltonians become, after identification of constants — 2&i = a, — \b 3 = /3, — 362 = 7, 5 2 = 
identical with the three-parameter Drach systems given in 20 . 

Let us now perform the transform between both families, i.e. between ([52]) and ([53]) . According to 
formula (|4ip in Theorem [7] the parameter independent parts of Hamiltonians transform as 

H = -W 5 -^R(F 7 )H = CH 

with R(F 7 ) = F' 1 and with 

_ ( vi 2) 



so that 



v 2 l 










-1 


) - ( -%i 


-0 


• 














4xj+xl J ■ 






XIX2 / 





c 

Moreover, the map (|43|) between solutions ([45]) of ([52]) and ([53]) becomes 

dt = -W^R-^F^dt = A T dt 

where 

_ 1. _ d/i _ / + ±a; 2 1 

G _ ^ -1*1*1 

and this map leaves the common level surface M a 5 invariant. Note also that now P(A, ot\, a 2 , 5i, 5 2 ) 
A _1 P(A, 5i, 5 2 , ai, a 2 ). The point transformation between fiat coordinates of both families is 

1 o 1 9 1 



xi = -2x 1/2 , x 2 = 2(y-x) 1/2 => x = -a; 2 , y = -x\ 



2 1 ^2 



2 " " 4 1 4 



a- 



14 



In our second example we will relate the HH family of separable potentials with a family of elliptic 
separable potentials. Both systems will belong to the same class of Stackel systems i.e. ji = Si for 
i = 1,2 now (so that W 7i a = —I and = Fg which is valid in any coordinate system) but they have 
different metrics. Let us thus first recollect some basic fact about generalized elliptic coordinates and a 
hierarchy of elliptic separable potentials [22] ■ Denote by (qi,...,q n ) the Euclidian coordinates on R™ and 
by (pi, ...,p n ) the conjugate momenta. The generalized Jacobi elliptic coordinates (Ai, A„) are defined 

1 + 4l t (,-^)- n ^ 1 (.-A')' (56) 
where /3, are nonzero different constants. Let us introduce the following abbreviations 

B(z) = f[(z-P j ), A(z) = f[(z-X j ), J ^L=B k (z) = -j2^P-^- i , (57) 

where are Viete polynomials with respect to its arguments. Then (|56[) takes the form 

n 

B{z) + l Y J B k {z)ql=k{z). (58) 
fe=l 

For z = A : 5(A) = 0, B fe (&) = <5 fcl B fc (/3 fc ), hence 



and then from (|5"T|) and 



Elliptic coordinates are separation coordinates for the following family of Benenti systems 

a(X) + ili A"" 1 + ... + H n = -i-B(A)^ 2 , 

where the family of elliptic separable potentials is given by cr(A) = A fc , fc e Z, and thus V( fc >(c?) = 
F fc (g)F(°), where F is given by ([291) and ([59]>. 

Let us relate two separation curves for n = 2, namely 

P(A, ft, a) = cr(A) + M + ft 2 - S" 1 (A)A(aiA + a 2 ) = -A/? (60) 

and 

P(X,h,a) = -B(X)a{X)X~ 1 - S(A)(5i + a 2 X^) + hiX + h 2 = -^B{X)fi 2 (61) 

where both families are now from Benenti class with (71,72) = {61,82) = (1,0) and with R(X) = 
— B{X)X~ X = —(A — Pi)(X — /?2)A _1 , while <r(A) is again assumed to be an arbitrary rational func- 
tion of A. For the extended Henon-Heiles system, when <j(A) = A 4 in ([60[) . in flat orthogonal coordinates 
from the previous example, an appropriate Hamiltonians are 



1 2 1 2 3 1 2 , 4(/? 1 +/3 2 )a: 2 -16/3i/3 2 x 1 

2»1 + 2 ?/2 + X 1 + - Xl X 2 + ^ - ^ _ ^ )(4/3 2 + Ahxi _ xl) < 



4x 2 2 + 



J? + 4&zi - x 2 )(4/? 2 2 + 4^ - x 2 )" 2 ' 



1 2 1 1 2 2 1 4^ + 

-xiy 2 gxayuft ^ -x 2 + _ ^ . 

4xixj + 4(01 + /3a)g| 

(4/? 2 + 4/3ia;i - x 2 )(4/3 2 + 4/3 2 X! - x 2 ) 



a 2 . 
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We pass now to the respective system from the family (IfTTl) . The transformation from Viete to fiat 
orthogonal coordinates is given by (|5§1) 

Pi (A) = -h -02 + \<ft + \ql, P2(X) = 0102 - \P2ql - \Piql 
The metric G and the Killing tensor K 2 of this system are 

1 \ ~ _( -p 2 + \q 2 2 -kiga 



a — — ^ 1 4^2 4' 

G ~ 1, 1 J' ^" -A + ig? 

so that the geodesic Hamiltonians of (fBTj) are 

^1 = \v\ + \v\, e 2 = ^(-/?2 + \q 2 2 )v\ + + \q 2 )pl - \q\qiV\V2- 



The matrix _F 7 becomes 



and hence 



7 I -0102 + \foq\ + \0iq\ 
T, ? 1 t7M ^ 2, 2\~ 02<ll+0iql 



where V^ CT ) = -B(F 7 )F^°\ so 

^ = + /? 2 3 9 2 2 ) + JtoA + /? 2 + /J 2 2 ) 9l <z 2 2 + ^Ctfrf + /? 2 2 <z 2 4 ) - ^(2A + /? 2 )<z 4 g 2 2 

- ^(A + 2/3 2 ) 9 2 g 2 4 - + fcq%) + ^(qUl + q 2 A) + ^qUI + ^(i s i+ 

V 2 (a) = j0i02(0Wl + Phi) - ^C 8 ? + + 2/?i/3 2 2 + PDqlql - ^0102{0A + foqt) 

4 16 16 

+ ^(01 + \0102 + PlMql + qlq 4 2)^0102(q! + ql) - ^(Wl + 02)q\q% 
According with (|48l) . we have now 

H = R{F 1 )H = —B(E y )F~ 1 H, 

where 

o/p s ( 01+02 + X! -l-Afc^X? 

1 7j ^-|4-A^ 2 l +fi 2 -^ihx l x-' 
while the reciprocal transformation (|49[) attains the form 

dt = R(F^)dt^ - ( || J * (62) 



where 



W ) = 



(fcql + 01<&)(±0102 ~ foq\ ~ 0iql)~ 1 ~[0hl + 0\ql){^2 ~ 2 q\ ~ Piql)' 1 

In (|62|) we used the inverse of formula (|49|) to avoid the complicated matrix R^ 1 (F 1 ). As before the trans- 
formation leaves the common level surface M a fi invariant. Also, this time we have P{\, a\,a 2 ,a\, a 2 ) = 
— B(A)A _1 P(A, ai, a 2 , ai, a 2 ). The point transformation between flat coordinates of both systems is 

Si - -01 -02 + \q\ + \q\, 4 = iq l + foq\ - i0x0 2 . 
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